Reconstructed Rough Growing Interfaces; 
Ridgeline Trapping of Domain Walls 
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We investigate whether surface reconstruction order exists in stationary growing states, at all 
length scales or only below a crossover length, Zrec- The later would be similar to surface roughness 
in growing crystal surfaces; below the equilibrium roughening temperature they evolve in a layer- 
by- layer mode within a crossover length scale Zr, but are always rough at large length scales. We 
investigate this issue in the context of KPZ type dynamics and a checker board type reconstruction, 
using the restricted solid-on-solid model with negative mono-atomic step energies. This is a topol- 
ogy where surface reconstruction order is compatible with surface roughness and where a so-called 
reconstructed rough phase exists in equilibrium. We find that during growth, reconstruction order is 
absent in the thermodynamic limit, but exists below a crossover length l^cc > ^r, and that this local 
order fluctuates critically. Domain walls become trapped at the ridge lines of the rough surface, and 
thus the reconstruction order fluctuations are slaved to the KPZ dynamics. 

PACS number(s): 64.60. Cn, 02.50.Ey, 05.40.-a, 68.35.Rh 




I. INTRODUCTION 



Equilibrium surface phase transitions have been a topic 
of research for several decades. Various types of critical 
behaviors are well established both in theoretical models 
and actual experiments. This includes surface roughen- 
surface melting |||,^, and surface reconstruc- 
5-0. Moreover, the competition between these 
phenomena leads to additional phases and phase tran- 
sitions, like disordered flat phases, preroughening tran- 
sitions, and reconstructed rough phases [p[-10|. Rough- 
ening induced deconstruction in Pt(llO) Mll| and pre- 
roughening induced deconstruction in Si(llO) type ge- 
ometries are other examples of this competition . 

The theory of dynamic non-equilibrium processes like 
surface growth has flourished during the last decade as 
well. Several new types of dynamic universality classes 
have been identified. KPZ type growth is one exam- 
ple |p^[-p2|. Unfortunately, in this area a gap seems 
to widen between theoretical and experimental interests. 
Theoretical oriented research tends to focus on universal 
aspects of these processes, such as the large scale prop- 
erties of growing surfaces in the stationary growing state 
and on how this state is approached in the asymptotic 
large time limit. Experimental oriented research tends 
to focus on more microscopic short distance aspects of 
growing surfaces, e.g., as encountered in actual epitaxial 
growth. 

One of the fundamental issues, relevant to both per- 
spectives, is whether any of the above equilibrium surface 
phase transitions persist in the stationary state of grow- 
ing interfaces. In this paper we address whether surface 
reconstruction order can exist during growth. 

This issue is related to the absence of surface rough- 
ening transitions in growing surfaces. Below the equi- 



librium roughening transition temperature Tr the grow- 
ing surface is rough at large length scales, but remains 
flat and grows layer-by-layer at distances shorter than a 
crossover length scale ^r, which varies with temperature 
and oversaturation. We review this briefly in section || 
in the context of elementary nucleation theory. 

Consider a surface that is flat and reconstructed in 
equilibrium at low temperatures. Below Tr it appears to 
grow within as flat in a layer-by-layer mode. More- 
over, below Tree (if Jlec < Jr) it appears as reconstructed 
if the new particles can find their proper reconstruction 
positions at times scales that are short compared to the 
rate at which a new layer is completed. Presume that 
this is indeed the case. The next, more intriguing ques- 
tion is whether ^rcc can be larger than Zr; i.e., whether 
rough growing surfaces be reconstructed? The compat- 
ibility of surface roughness with surface reconstruction 
was addressed in the context of equilibrium phase transi- 
tions several years ago. The answer depends on intricate 
details of the surface topology. For example, in missing 
row reconstructed (MRR) (110) facets in FCC crystals, 
like Au and Pt, roughness is incompatible with recon- 
struction order, and the surface roughening transition 
must destroy the reconstruction simultaneously In 
such geometries, the reconstruction order can not exist 
in growing surfaces beyond the roughness length scale 
either, and ?rec < Ik- 

Surface roughness and reconstruction are compatible 
with each other in other crystal structures. Simple cubic 
(SC) MR reconstructed (110) facets are an example. In 
equilibrium, they can roughen before the reconstruction 
order deconstructs, Tr < Tiec- The intermediate phase is 
known as a reconstructed rough phase 0,|l^. For those 
surfaces it might be possible to observe genuine decon- 
struction type phase transitions in growing surfaces. Or, 
if not, the surface reconstruction can at least persist well 
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beyond the roughness crossover length scale, Zrcc > ^r, 
and will be limited by an independent mechanism. These 
issues are the topic of our research reported here. 

In section |l| we review rough versus layer-by-layer 
growth in surfaces, and in section III the basic properties 



of equilibrium reconstructed rough phases. Next, in sec- 
tion [V , we start to focus on the reconstruction versus dy- 
namic roughness issue, and then, in section ^ we choose 
a specific type of reconstruction and a specific type of sur- 
face growth dynamics to study it quantitatively by means 
of Monte Carlo (MC) simulations. The model must be as 
simple as possible, avoiding secondary effects that might 
obscure the central issue. Our choice is the so-called re- 
stricted solid-on-solid (RSOS) model with negative step 
energies, which describes a simple cubic checker board 
type reconstruction, and KPZ type growth. The MC sim- 
ulation re sults are presented in section and analyzed 
in section VII . Finally, in section VIII| wc summarize our 
results. 



II. ROUGHNESS IN GROWING SURFACES 



The topic of this paper is whether surface reconstruc- 
tion order can exist during growth, but as a start it is use- 
ful to review briefly the related issue of dynamic surface 
roughness from a long and short length scale perspective. 
Elementary nucleation theory suffices for this purpose. 
Equilibrium crystal surfaces undergo well defined rough- 
ening transitions from macroscopic flat to macroscopic 
rough. On the other hand, growing surfaces are theoreti- 
cally "always rough" |l^,|lj]. This seems at odds with 
practical reality, where surfaces appear to grow quite 
differently below and above the equilibrium roughening 
temperature Tr. Above Tr they are rough (dynamic 
roughness) while below Tr they seem flat (layer-by-layer 
step- flow growth). Above Tr the growth velocity Vg is 
proportional to the oversaturation Vg ~ Afi, while below 
Tr it is inversely proportional to a nucleation time scale 
Vg ^ with exp {-arj'^ / {AfiksT)) |jl||. 77 is the 

equilibrium step free energy. As a result, crystal growth 
shapes have sharp angles, in which many facets, including 
all that are above their Tr, are missing. This apparent 
difference in growth mechanism is one of the most use- 
ful experimental tools to locate equilibrium roughening 
transitions in crystal facets. 

The origin of the exponential factor in t is the ex- 
istence of a nucleation barrier for creating a terrace of 
height h ^ h + 1 below Tr . The edge (step) free en- 
ergy loss term (proportional to 77 times the circumfer- 
ence) competes with the surface energy gain term (pro- 
portional to A/Lt times the terrace area). The nucleation 
barrier vanishes when the step free energy 77 vanishes, 
i.e., at Tr. After a new terrace larger than the nucle- 
ation thresholds is nucleated with an exponential small 
probability, it spreads out fast by particle adhesion at 
its edge into a macroscopic domain, until it merges with 



other spreading terraces that have nucleated in the mean 
time, and thus complete the new surface layer. However, 
new terraces are nucleated on top of spreading terraces 
as well. This nesting effect, together with the spatial 
fluctuations of nucleation events leads to a loss of a well 
defined (length scale free) global reference surface level. 
This means that although at small enough length scales 
the surface looks flat and seems to grow layer-by-layer, 
at large length scales it is rough. 

There is no phase transition between the layer-by-layer 
and rough growth regimes, only a characteristic crossover 
length scale. The latter is of order ^r — VgT, with Vs the 
step velocity (determined by the particle deposition rate 
at the step edge) and r the above time scale at which ter- 
race nuclei are being created. Surface flatness cannot be 
maintained during growth over large length scales, but 
at small oversaturations (A/i) and sufficiently below Tr 
(large step free energies 77) the growing surface can ap- 
pear to be flat for all practical purposes, over any typical 
experimental length scale. 

The same type of issues arise in our study concern- 
ing the compatibility of surface reconstruction order with 
growth dynamics. First we address whether surface re- 
construction order can persist during growth at macro- 
scopic length scales (the thermodynamic limit); and, if 
not, whether it might still exist in a practical sense within 
a characteristic length scale, ^rcc below the equilibrium 
reconstruction temperature T^-cc- 



III. RECONSTRUCTED ROUGH EQUILIBRIUM 
PHASES 

Surface reconstruction is conventionally associated 
with flat interfaces. However, surface roughness not 
necessarily destroys the reconstruction order. A rough 
but still reconstructed surface is in a so-called recon- 
structed rough (RR) phase. The equilibrium versions of 
RR phases were studied theoretically some years ago in 
the context of the competition between surface roughen- 
ing and reconstruction in MR reconstructed FCC (110) 
facets [^,0. The topological details of those FCC sur- 
face prevent the existence of RR phases, implying that 
in Pt(llO) the surface roughness and deconstructs simul- 
taneously 10] as observed experimentally in Pt(llO) [ pT| . 
This implies immediately that during growth reconstruc- 
tion order is limited to the roughness crossover length 
scale, ^rcc £ ^R- The same theoretical studies also identi- 
fied other surface geometries where RR phases do exist. 
For those /roc is not limited by l^. In this section we 
review the basic properties of RR phases, using as ex- 
amples checkerboard and MR type reconstructed simple 
cubic stackings. 

To avoid confusion, it is useful to distinguish between 
misplacement and displacement type reconstruction [ p^ . 
In misplacement reconstructions, particles have moved to 
different solid-on-solid type stacking positions, or are re- 
moved altogether, compared to the unreconstructed flat 
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surface structure. The checker board reconstruction in 
Fig.|l| and also the more realistic MR type reconstruc- 
tions are examples of this. The average surface height 
has changed by half a unit, /i — > /i — i. In displacement 
reconstructions the atomic stacking does not change. In- 
stead, the atoms are merely elastically distorted at the 
surface with a commensurate or incommensurate period 
compared to the bulk. Misplacement type reconstruc- 
tions are more likely to disorder at temperatures near Tr 
than displacement type reconstructions. For clarity we 
focuses here on misplacement reconstructions. 

The definition of the reconstruction order parameter is 
at the core of RR phases [0,0. In checker board and 
MR reconstructed SC (110) facets , the reconstruction 
order can be formulated in two distinct ways. One for- 
mulation keeps track of whether the black or white fields 
(even or odd rows) are on top. The other measures it 
in terms of anti-ferromagnetic order in the parity type 
Ising variables Sr — eKp^inhj.), with hj. = 0, ±1, ±2, • • • 
the surface height at site r, see Fig.^j. These two formu- 
lations might seem equivalent in fiat surfaces, but they 
are not in the presence of roughness. 

The compatibility of surface reconstruction with sur- 
face roughness depends on topological properties of step 
and domain wall excitations; on how they affect the two 
versions of the order parameter. Fig.^ shows in cartoon 
style a cross section of the reconstructed surface, and also 
domain wall and step excitations. The domain wall in (b) 
does not change the surface height. Notice that both or- 
der parameters change sign. Across the step in (c), from 
left to right, the even-odd order parameter changes sign, 
but the parity order is unaffected. At the step in (d) the 
opposite happens. 

These two types of steps are the only topologically dis- 
tinct ones that are possible; (c) couples only to the even- 
odd row type order parameter and (d) only to to the 
parity version. 

It is possible to construct many more step and domain 
wall structures that look locally different from the ones 
in the figure, but those induce the same change in height 
and/or reconstruction order(s) and therefore are from a 
topological point of view identical to the ones in the fig- 
ure. Notice also that the excitations in (b)-(d) are related 
to each other in the sense that any of the three can be 
interpreted as a bound state the two others. Elastic sur- 
face deformations in the actual atomic positions near the 
surface and additional ones near the steps and domain 
walls, influence the local internal structure of steps and 
domain walls, but do not affect those topological fea- 
tures, and therefore need not explicitly be represented 
in the following discussion. (They certainly renormalize 
the step and domain wall energies and the interactions 
between such surface excitations.) 

The fate of reconstruction versus roughness depends 
on the energies of these steps and domain walls, includ- 
ing the kink energies. They set the scale of the meander 
type entropy and therefore the temperature dependence 
of the step free energies. If the domain wall free energy 



vanishes first, the surface remains fiat but the reconstruc- 
tion vanishes, Tree < Tr- In case the free energy of one of 
the two types of steps vanishes first, the surface enters a 
reconstructed rough phase, Tr < Ticc. At the roughening 
transition one of the two reconstruction order parameters 
vanishes, but the other type of order remains. So there 
exist two topologically distinct types of RR phases. (No- 
tice that only the one with the parity type order is readily 
observable by, e.g., conventional X-ray diffraction.) 

In the RSOS model below, the RR phase has parity 
order, i.e., the step free energy of the (c) type steps is 
zero, but walls and (d) type steps still have non-zero free 
energy. We will refer to those excitations as "loops of 
zero's", because in the rough surface they show up as 
contours across which the height change is zero, dh = 0. 
The deconstruction transition (inside the rough phase) 
takes place at the temperature where the surface tension 
of the loops vanishes. In equilibrium that turns out to 
be an ordinary Ising transition. This conclude our brief 
review. For more details we refer to Refs. and |l0[| . 



IV. RECONSTRUCTED ROUGH GROWTH 

Let's focus now on surface growth. Only in surfaces 
where equilibrium RR phases are topologically possible, 
can the surface reconstruction length scale Irec exceed the 
onset of dynamic roughness length scale Zr. Moreover, 
it is quite possible that the reconstruction order persists 
over all length scales (l^cc ^ oo), such that a genuine 
dynamic deconstruction phase transition takes place in 
the stationary state of the growing surface, just like in 
equilibrium. 

For comparison, imagine a two dimensional (2D) lat- 
tice with on each site an height variable and an Ising 
spin degree of freedom (representing the reconstruction 
order). This leads to two coupled master equations, one 
for surface growth, e.g., KPZ type dynamics, and another 
for the reconstruction order, e.g., Glauber type Ising dy- 
namics. In equilibrium surfaces, the coupling between 
the two sectors is weak, to the extend that the recon- 
struction transition in the Ising sector and the roughen- 
ing transition in the height variable sector do not interfere 
with each other 0j2^. The central issue is whether and 
how this coupling changes during growth. The Ising dy- 
namics itself, is blind to the growth bias. If the coupling 
between the two sectors is remains weak, the Ising spins 
can still reach the Gibbs equilibrium state and undergo 
a conventional equilibrium reconstruction transition. 

Coupled master equations of this type have been stud- 
ied recently in the context of specific ID growth mod- 
els. Those display strong coupling between the Ising 
and roughness degrees of freedom, such as growth be- 
ing pinned down by Ising domain walls p^-|2^]. Pinning 
favors to spontaneous facetting. In our 2D model, we 
observe different effects, besides the obvious fact that in 
ID equilibrium reconstruction order can not exist. 
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V. RESTRICTED SOLID-ON-SOLID MODEL 

The 2D restricted solid on solid (RSOS) is one of the 
work horses of surface physics research. Integer valued 
height variables hr = 0,±1,±2,-- - are assigned to a 
square lattice and nearest neighbor heights are restricted 
to differ by at most one unit, dh =^ 0, ±1. The energy 



^ = E 



(1) 



(r.r') 



depends only on nearest neighbor interactions. We use 
dimensionless units, K = J/k-QT. The if > side 
of the phase diagram contains a conventional equilib- 
rium surface roughening transition p3t . Moreover, the 
non-equilibrium version has been studied extensively for 
A' > as well, because it is a natural lattice realization 
of KPZ growth 

For if < 0, the model contains one of the simplest ex- 
amples of an equilibrium RR phase , and is probably 
the most compact formulation of the coupling between 
Ising and surface degrees of freedom. The c?/i = ±1 steps 
are more favorable than flat dh = segments. At zero 
temperature, K — > — cxd, the dh = states are frozen 
out, and the model reduces to the so-called body cen- 
tered solid on solid (BCSOS) model, but in this version 
it lacks step energies, which means that the surface is 
rough even at zero temperature. The surface is rough, 
but since nearest neighbor heights must differ by one, all 
heights on one sublattice are even and odd at the other, 
or the other way around. This two-fold degeneracy rep- 
resents the checker board type RR order. The staggered 
magnetization, defined in terms of the parity spin type 
variables Si = exp^inhr), is non-zero. 

The dh = excitations that appear at T > form 
closed loops and behave like Ising type domain walls. 
The reconstruction order changes sign across such loops. 
Their sizes diverge at the equilibrium deconstruction 
transition Kc — —0.9630 |23|. (Determined by trans- 
fer matrix finite size scaling techniques). The Ising and 
roughness variables couple only weakly. Numerically, all 
reconstruction aspects of the transition follow conven- 
tional Ising critical exponents. Moreover, the thermo- 
dynamic singularities in the Ising sector affect only the 
temperature dependence of the surface roughness param- 
eter Kg, defined in terms of the height-height correlator. 



{{hr+r„~hr,f)^{7TKG)-'Hr). 



(2) 



The continuum limit analysis confirms these numerical 
results. The point in the generalized phase diagram 
where the Gaussian (height) and Ising degrees of freedom 
decouple is a stable renormalization type fixed point . 

We study this same RSOS model in the presence of a 
KPZ type growth bias. In the MC simulation, we first 
select an update column and next whether a particle de- 
position or evaporation event will be attempted. The 
move is rejected if it would result in a violation of the 



RSOS condition, dh = 0,±1. If allowed, it will take 
place with probability P — mm{p,pe~^^') in case of de- 
position, and with probability P = imn{q, qe^^^^ ) for 
evaporation. Without loss of generality we can choose 
p + q — 1. At infinite temperature {K = 0) and depo- 
sition only {q ~ 0) the model reduces to the well known 
Kim-Kosterlitz pQ] model for KPZ type growth. 

We will present only our MC results far from equilib- 
rium, i.e., at g = with deposition only. We observe no 
qualitative differences closer to equilibrium, < p < 1, 
but the interpretation of the data becomes increasingly 
obscured (as expected) by (conventional) crossover scal- 
ing from the equilibrium deconstruction phase transition. 

At low temperatures, K — » — oo, the Metropolis dy- 
namics slows down considerably. The rejection rate be- 
comes high and the density of active sites becomes low. 
Therefore we employ the following rejection free algo- 
rithm. During the MC simulation we keep a list of ac- 
tive sites, i.e., sites where particles can deposit with- 
out violating the RSOS condition. They are grouped 
in j = 1, ■ • ■ , 5 sets, according to the five distinct en- 
ergy changes A.Ej that can occur during deposition. 
First we preselect one of those 5 sets, with probability 
{PjNj)/{^j PjNj), where pj = min(l, e"'^-^^) and Nj is 
the number of sites of type j. Next, a particle is ran- 
domly deposited at one of the sites in that specific set 
j. Rejection free procedures like this upset the flow of 
time. To restore proper time, we increase the MC time 
during each update step by 1/p x 1/Nj. We checked 
explicitly that this reproduces the correct value for the 
KPZ dynamic exponent z — 8/5 pO|-p^ a,t K ^ 0; we 
find z ~ 1.6 ±0.1. 

The above algorithm resolves the slowing down prob- 
lem in the actual MC simulation, but does not address 
its origin. In the limit K — > —oo the RSOS model re- 
duces to the BCSOS model, with d/i = ±1 at all bonds. 
The dh = loops are frozen out completely. In BCSOS 
type KPZ growth dynamics, 2 particles are deposited 
at ones in the form of vertically oriented bricks, oth- 
erwise a "forbidden" configuration with dh = would 
arise. In the K < RSOS model at very low tem- 
peratures the same event is achieved as a 2-step two- 
particle process, by the deposition of a second particle 
at the same site soon after the first one. The prob- 
ability for deposition of the first particle is equal to 
p = exp{2K). The second particle deposition on top 
of it happens with probability p — L^'^). This implies 
that the time clock in the RSOS model runs slower by a 
factor r = exp(2i(:)(l -I- 4exp(i4:) -I- ...). 

A final remark about surface roughness. In normal sur- 
faces, the equilibrium roughness increases with tempera- 
ture; due to the fact that meander type entropy renormal- 
izes the step energy into a reduced step free energy Q . 
In our model, surface roughness evolves the opposite way; 
it decreases with temperature. The surface is less rough 
at infinite temperature K = Q than in the zero tempera- 
ture limit K — > — oo. A high temperature RSOS surface, 
with dh = 0, ±1 is obviously less rough than a BCSOS 
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surface, with only dh = ±1. Recall that this BCSOS 
model lacks step energies, such that it is just as rough 
at T = as at T ^ oo. From the BCSOS perspective 
the thermally excited dh = loops stiffen the surface, 
and give rise to an inverted roughness versus tempera- 
ture profile. On the one hand, this is an interesting phe- 
nomena in its own right. Moreover we could fine tune it 
by introducing next-nearest neighbor interactions, since 
they represent BCSOS type step energies. On the other 
hand, this effect is unlikely to affect the central question 
we want to address (how do roughness and reconstruction 
degrees couple during growth) and therefore we choose 
not do so in this study. 

VI. RECONSTRUCTION DURING GROWTH 

We search for reconstruction order as function of tem- 
perature, for —oo < K < 0. The susceptibility type 
parameter p7| ], 

X^L\{m')^{\m\f) (3) 
of the reconstruction order parameter, 

^ ^ (^(^^ly+y g^^h{x,y)'^ (4) 

is shown in Fig.|^ for the stationary state of the grow- 
ing surface, as function of K for different system sizes 
L^. The sharp maxima seem to confirm the existence 
of a dynamic surface reconstruction transition into a RR 
phase. However, several features are very different from 
equilibrium. The peak height diverges as ^ ^ . i.e., 
stronger than at the equilibrium transition point where 
it scales as x ~ L'^^" . This could be a signal of a first 
order phase transition. However, the peak position does 
not converge to a specific critical point Kc- Instead it 
keeps shifting with lattice size. It scales logarithmically, 
as ii:poak(i) =i -A\n{L/Lfi) with A = 0.77 ± 0.05 and 
Lq = 2.2 ±0.2. 

Next, we monitor in detail the reconstruction order 
parameter m near and below the equilibrium Kc as func- 
tion of time. It behaves similar as in conventional spon- 
taneously ordered phases, but flip-flops more frequently 
than justifiable from finite size effects alone. Moreover, 
the fluctuations in m within each phase are too strong. 
Fig.^ quantifies this in terms of a histogram of the num- 
ber of times a specific value of m appears in a typical time 
series. The distribution has two distinct peaks, suggest- 
ing the presence of spontaneously broken reconstruction 
order, but the tails have a power law shape instead of 
the exponential form mandatory for a spontaneous bro- 
ken symmetry. 

Power laws are the hallmark of critical fluctuations. 
So, quite surprisingly, it appears as if the RR order is 
critical at low temperatures for all K < Kp^ak- Instead 
of an isolated critical point, we seem to be dealing with 
a critical phase. 



VII. LOOPS TRAPPED ON RIDGE LINES 



The surprising critical fluctuations in the reconstruc- 
tion order parameter can be traced to the following loop 
dynamics. Consider a typical configuration at very low 
temperatures. Fig .|| shows an example ||2^. The surface 
is in an almost pure BCSOS type dynamic rough station- 
ary state (with d/i ± 1), and contains only a few dh — 
loops separating surface areas of opposite checker board 
type RR order. 

The typical life cycle of such a loop runs as follows. 
It is nucleated in a valley bottom. Next it runs up hill, 
grows in diameter, encompassing the entire valley, until 
it becomes trapped on a ridge line. There it lingers un- 
til another loop annihilates it, or when the KPZ surface 
fluctuations to which it is slaved shrink it back to zero. 

Fig.^(a) represents a cross-section of the 2D rough sur- 
face near a valley. It shows a domain of opposite re- 
construction inside an otherwise perfectly reconstructed 
rough configuration. The two flat segments are the lo- 
cations where the domain wall loop intersects the cross 
section. In equilibrium, the loop fluctuates with equal 
probability up and down the slope because depositions 
and evaporations are equally likely. A growth bias breaks 
this symmetry, the loops move more likely upwards than 
downwards, see Fig.^(a). This upward drift is the driving 
force responsible for the trapping of loops at ridge lines, 
and thus creates a strong coupling between the roughness 
and reconstruction degrees of freedom, unlike equilibrium 
where they effectively decouple. 

A few comments on the topology of ridge lines in rough 
surfaces might be useful. Imagine a rolling ball in this 
landscape, like in the well known analogy with renormal- 
ization flow in statistical physics. Presume strong friction 
such that the velocity is proportional to the force, i.e., the 
gradient of the slope, at all times. The hill tops are the 
completely unstable "fixed points". The valleys are the 
attractors. The ridge lines form the water sheds between 
valleys. Every ridge line runs from an hill top to a saddle 
point. From each hill top an arbitrary number of ridge 
lines can emerge, but only two ridge lines can end at each 
saddle point (at opposite sides of the single direction in 
which the saddle point attracts). So the ridge lines form 
a network, and since none of them can not stop in midair 
it is a closed network. The KPZ rough surface is scale 
invariant, which means that this ridge line network has 
fractal properties. 

Ignore for the time being the scale invariant aspects 
of the network. Imagine a landscape consisting of deep 
smooth valleys surrounded by ridge lines; unlike the real 
rough surface where every deep valley consists of col- 
lections of sub valleys. The life cycle of a macroscopic 
loop in this surface starts with the nucleation of a new 
seedling-loop at the floor of the valley and its rise along 
the slopes, during which it grows into a macroscopic ob- 
ject. The only loops of interest are those nucleated at 
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the valley bottom and then run up-hill encompassing the 
entire valley. Only those loops are topologically trapped 
and stable. Loops nucleated on the slopes annihilate by 
stochastic fluctuations before becoming macroscopically 
large. The same is true for loops nucleated out of the 
valley bottom but running up-hill on one slope segment 
only. 

The rise of a seedling loop out of the valley bottom into 
a macroscopic object, is a very fast process. Almost no 
MC moves that make the loop grow and rise are rejected; 
energy barriers are rare, because the length of the loop 
(its energy) increases uniformly. Compared to this, the 
nucleation frequency in the valley bottom is very small. 
This means that the time scale at which a macroscopic 
loop emerges out of the valley is limited by the nucleation 
time scale Tn and independent of the valley size. 

To measure t,i we prepared a surface in the BCSOS 
KPZ stationary state and measure (at a very low tem- 
perature, K <^ Kc) the intervals between macroscopic 
loop events. Numerically we find t„ ~ exp(— aX) (mea- 
sured in BCSOS time units) with a = 3.0 ± 0.1. 

This agrees qualitatively with the following estimate. 
The deposition of the first particle in the valley bottom 
occurs with probability p = L~^e^^ . This creates a 
fledgling loop, but one that is indistinguishable from the 
intermediate state in an elementary BCSOS type growth 
event (where a second particle is dropped on top of it 
with probability p = L^'^). The loop grows when the 
next particle is dropped not on top but next to the pre- 
vious one. That happens with probability p = L~^e^ . 
The nucleation threshold diameter Ic is reached when the 
loop growth and BCSOS growth become distinguishable, 
i.e., when the annihilation of a loop requires the creation 
of a new well distinguishable loop inside it. That happens 
at about ~ 7, see Fig.^. The time scale at which that 
stage is reached is approximately t ~ L^^e""^^ (in BC- 
SOS time units), which is of the same order of magnitude 
as the above numerical nucleation time scale. 

The loop rises out of the valley until it becomes 
trapped on the ridge line that borders this valley to ad- 
jacent ones. From there on the loop is slaved to the 
growth fluctuations of the surface. Valleys grow and 
shrink (without bias), open- up, fill- up and merge. The 
loop has to follow this dance of the ridge line until a 
new loop nucleates out of the valley and annihilates it, 
or when the encircled terrain happens to shrink to zero 
(fills- up) by surface growth fluctuations. 

We expect that the life time Tz{L) of a ridge line of size 
L in a growing surface, scales as a power law, Tz ^ L^, 
with z the dynamic exponent of the surface roughness 
degrees of freedom (KPZ like in our model). To test this, 
we measure the decay times of large macroscopic defect 
loops (of about half the lattice size) as function of L, at 
low temperatures K <C Kc- The data in Fig.|^ collapse 
indeed onto one universal curve after a rescaling of time 
by Tz « L^. The collapse fits best at z = 1.7 ± 0.1 (in 
BCSOS time units), which is consistent with the known 
KPZ dynamic exponent z = 8/5 |^-|2|. 



The ridge line fluctuations are responsible for the 
power law tails in the time distribution of RR order, 
Fig.^. Those critical fluctuations only show up below a 
characteristic length scale Zioc, where the nucleation time 
scale Tn ~ exp(— aif) is larger than the surface growth 
time scale, L^. A simple estimate for l^^c follows 

from equating the two time scales, l^ec exp(jK). 

The peaks in the susceptibility, in Fig.||, reflect this 
crossover length l^cc- Recall that the peak shifts loga- 
rithmically. By setting t„ = Tz we obtain the same log- 
arithmic behavior, Kc = —■^\n{L/Lo). The prefactor is 
too small by about 30%, but this is not a surprise because 
the estimate is rather simple minded. It ignores for ex- 
ample the self similarity of the rough surface. Consider 
a sub valley adjacent to an already trapped loop. Sup- 
pose a new loop nucleates out of this subvallcy. The loop 
segments annihilate each other in pairs. The net effect of 
this nucelation event is therefore that the trapped loop 
jumps across the sub valley. It now follows the comple- 
mentary segment of the ridge line that encircles the sub 
valley. Such events renormalize t^, in particular near Kc- 



VIII. CONCLUSIONS 

In this paper we study the compatibility of surface re- 
construction and surface roughness during growth. There 
are several possibilities. 

In surfaces where reconstructed rough (RR) phases are 
topologically forbidden, like missing row reconstructed 
FCC(llO) facets, the reconstruction order can not exist 
on a global scale in the stationary growing state. It can 
appear only locally within the crossover roughness length 
scale /r within which the surface grows in a layer-by-layer 
fashion, i.e., Ij-cc < ^r- 

The reconstruction length scale Zroc can exceed Zr, only 
in surfaces where equilibrium reconstructed rough phases 
are topologically possible, and those surfaces could in 
principle even display genuine deconstruction type phase 
transitions in the stationary growing state. 

We address this issue in the context of KPZ type dy- 
namics, in the RSOS model with negative coupling con- 
stant K < 0, which in equilibrium has a checker board 
type RR phase and a true deconstruction phase transi- 
tion inside the rough phase. We find that the stationary 
growing rough state lacks true macroscopic RR order; 
Ircc remains finite. Moreover, we identify the mechanism 
that sets the temperature dependence of Zioc- 

The fundamental features are an upward drift of the 
reconstruction domain wall loops and their trapping at 
the ridge lines of the surface. There, the loops are slaved 
to fluctuations of the surface growth dynamics, /^ec is set 
by the competition between two life times: the nucleation 
time scale of a new loop out of the valleys (annihilating 
existing trapped loops) and the time scale tkpz ^ L'^ 
at which a ridge line of radius L vanishes due to surface 
growth fluctuations. 
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At length scales smaller than l^cc ^ exp{jK), the sur- 
face appears as reconstructed rough, and the life time 
of the loops is determined by the KPZ growth dynamic 
fluctuations. The latter follow power laws. This man- 
ifests itself in critical fluctuations in the reconstruction 
order at length scales smaller than ^rcc- In x-ray diffrac- 
tion from such a growing interface, one would observe 
not only power law shaped peaks associated with the 
surface roughness, but also at temperatures where l^cc is 
larger than the coherence length of the surface, power 
law shaped reconstruction diffraction peaks. 

At length scales larger than ?rec, the surface appears as 
unreconstructed rough. Loops of that size die by nucle- 
ation of new loops instead of KPZ surface fluctuations, 
and they are not trapped anymore, because loop seg- 
ments can hop across sub valleys of size I > l^^c by means 
of nucleation of new loops in sub valleys. 

In our study we chose to focus on KPZ type surface 
growth dynamics, but we have good reasons to expect 
that the trapping of domain walls on ridge lines is a 
common phenomenon. In general, the quasi-critical fluc- 
tuations will reflect the dynamic exponent of whatever 
growth dynamics is applicable. In recent studies of ID 
models with KPZ and Ising type coupled degrees of free- 
dom the Ising defects become trapped in valleys and 
canyons and thus pin-down the growth 25|. We ex- 
pect that a tendency towards facetting instead of ridge 
line trapping can also be realized in our 2D model by 
varying the local growth rates. This research is supported 
by the National Science Foundation under grant DMR- 
9985806. 
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FIG. 1. Checker board type misplacement surface recon- 
struction 




FIG. 2. (a) A perfect reconstructed surface, (b) A do- 
main wall; both order parameters change sign, (c) A step 
where only the even-odd row order changes sign, (d) A step 
where only the parity order changes sign. 
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FIG. 4. Histogram (insert) of the reconstruction order pa- 
rameter, m, at L = 32 and K — —3.2 from 2^** data points 
using AM = 0.01 as bin width. The tails about the peaks 
at m = ±1 scale as power laws (main frame) with exponent 
-0.9 ±0.1. 
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FIG. 3. Reconstruction order susceptibility x ^ function 
of temperature at system sizes L — 8-64. The data collapses 
onto a single curve by the shift K' = K — 7i'peak(i), with 
ifpoak(L) = -0.771n(L/2.2). 



FIG. 5. A typical low temperature configuration of the 
growing surface with one large loop trapped at a ridge-line. 




FIG. 6. (a) One dimensional cross section of the surface 
near a valley with two loop segments. On the slope, a (d) 
are the only active adsorption (desorption) sites. The domain 
walls always move upwards during adsorption, (b) A loop of 
size of Ic nucleated at the bottom of a local valley. Gray and 
white sites have different surface reconstruction parity order. 
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FIG. 7. Histogram of the decay time of a trapped loop at 
K = —6.0. The data collapses by rescaling time by a factor 
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